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Abstract 

We investigate the dynamics of interacting magnetic flux lines driven 
by an external current in the presence of linear pinning centers, arranged 
either in a periodic square lattice or placed randomly in space, by means 
of three-dimensional Monte Carlo simulations. Compared to the non- 
interacting case, the repulsive forces between the vortices reduce the criti- 
cal current J c , as determined from the depinning threshold in the current- 
voltage (I-V) characteristics. Near the depinning current J c , the voltage 
power spectrum S(u)) reveals broad-band noise, characterized by a l/ui a 
power law decay with a < 2. At larger currents the flux lines move with an 
average velocity v cm . For a periodic arrangement of columnar pins with a 
lattice constant d and just above J c , distinct peaks appear in the voltage 
noise spectrum S(lj) at u> ~ v cm /d which we interpret as the signature of 
stick-slip flux line motion. 



1 Introduction 



Magnetic vortices in the mixed state of high-T c superconductors in the presence 
of defects — either localized point pinning centers, such as oxygen vacancies 
or extended columnar pins introduced into the material via high-energy ion ir- 
radiation — offer a complex system which has attracted the attention of both 
experimentalist and theorists from various disciplines of physics in recent years 
1 . This area of research is also inspired by important technological questions 
triggered by the search for effective pinning centers in order to minimize dis- 
sipation when current is passed through the superconductor 1 . From a more 
fundamental viewpoint, one is interested in the non-equilibrium dynamics of 
a disordered flux line system. When subject to an external electric current, 
non-trivial steady states are generated which can be directly probed experi- 
mentally. For example, the current-voltage (I-V) characteristics pQ, the voltage 
noise power spectrum 0II3E], or the local flux density noise [3] measured in 
experiments bear important signatures of the different non-equilibrium steady 
states. At low currents, transport can take place due to thermally activated 
'tunneling' of flux lines through the defects £Q, while at larger currents various 
intriguing steady states have been proposed, such as plastic flow where flux lines 
move in channels [o] . With increasing current the flux lines tend to reorder and 
move in a lattice (moving Bragg glass for point vortices [HE]), or the channels 
in the plastic flow regime can organize in a periodic structure (smectic phase) 
in the transverse direction of the external current \J\ . 

The presence of such steady states is reflected in the voltage noise spectrum; 
e.g., reordering would give rise to characteristic 'washboard' noise features. Re- 
cent experiments have in fact demonstrated the presence of both narrow- and 
broad-band noise in the measured voltage and local flux density spectra 01310]. 
Our goal is in addition to study how these noise spectra depend on the type of 
defects (correlated or point-like) present in the sample. However, since analyt- 
ical studies of the vortex dynamics are usually limited to asymptotic regimes 
IB] E] where the current J is cither much smaller or much larger than the 
critical current J c , one needs to resort to either Langevin molecular dynamics 
simulations 1101 1111 IT2"] or Monte Carlo techniques ^3 E] to study the vari- 
ous non-equibrium steady states which occur in the intermediate current regime. 
Owing to their complexity, these numerical simulations are largely restricted to 
two dimensions |51 ll0llTTlll2j and very few three-dimensional studies have been 
reported [131 1141 IT51 1161 117| . In this article we present the results of fully three- 
dimensional Monte Carlo simulations of interacting magnetic vortex lines in the 
presence of columnar defects. We investigate the effects of the repulsive forces 
between the flux lines and their interactions with extended defects in the I-V 
characteristics and the voltage noise spectrum, focusing on those aspects which 
are inaccessible by two-dimensional simulations. 

This article is organized as follows. The model is introduced in Sec. II and 
the Monte Carlo scheme is described in Sec. III. In Sec IV we present the results 
of our simulations for the I-V characteristics and voltage noise power spectra. 
Some conclusions are offered in Sec. V. 
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2 The model 



We consider a system consisting of n flux lines in the London approximation in 
the presence of N p defects. Our description utilizes an effective Hamiltonian 6 

+ ± e int [ L dsKo( {rixis) r rj±is)l ) +J2 f L dsV^is)] , 

i,j = 1;i^£j ' i=l ' 

where i"j(s) = [rjj_(s), Zi(s)] describes the configuration of the ith flux line in 
three dimensions, and the crystalline axis c of the superconducting material 
(as well as the external magnetic field) is oriented parallel to z, in a sample 
of thickness L. The first term gives the line tension energy of a magnetic flux 
line. This linear elastic form of energy holds good as long as \dr±(z)/dz\ < I/F 
[T%] , where F denotes the anisotropy ratio ^ |Hj ■ The line stiffness is given by 
e~£o ln(A b/Coh) DIE], wrtn £o = (0o/47rA a b) 2 (4>o = hc/2e is the magnetic flux 
quantum). A a fc and £ a f, denote the penetration depth and the superconducting 
coherence length in the ab plane, respectively. For high-T c materials, r 3> 1 pQ, 
giving the quadratic form of the line tension energy a wider range of validity. The 
repulsive interaction energy between the flux lines is approximated to be local 
in the z coordinate, in accord with the London limit, ft is represented by the 
second term in (Q, where e int = <\>\ / '(8n 2 X 2 ab ) and K (x) is the modified Bessel 
function which varies as — In a: when i<1 and ~ x~ 1 l 2 e~ x as x 3> 1. We model 
the pinning potential as a sum of N p independent potential wells, Vp[r(s)] = 

J2k=i UQ(r p — |r(s) — |), where r p is the radius, 8(x) denotes the Heaviside 

(v) 

step function, and the indicate the N p locations of the pinning sites. The 
free energy Fl(T) per unit length is then defined through exp[— /3Fl(T)L] = 
J Yli—i D[ri] exp (— f3H[ri\), where (3 = (ksT)^ 1 . We study the dynamics of 
the flux lines when an external current J is applied through the system which 
produces a Lorentz force = (l/c)z x J = (J/c) x, per unit length of each flux 
line. We resort to a three-dimensional Monte Carlo (MC) simulation to deal 
with the system in an intermediate range of driving currents where the system 
appears to be intractable by means of analytical methods. 



(1) 



3 Monte Carlo simulation 

In our Monte Carlo simulation each flux line r(s) is modeled by N = L/ao points 
where the ith point is located at r(i) = [x(i),y(i), z(i)] and interacts with its 
nearest neighbors via a simple harmonic potential | 2(j)=i-i i+i [F — 
rj_((j))) 2 + (\z(i) — z((j))\ — a ) 2 ]. The line is placed in a box of size L x xL y x L z 
with periodic boundary conditions in all directions. The interaction energy 
between a pair (I, m) of flux lines is given by e int J2i=i K (\r±i(i)-r± m (i)\/ X ab ). 
To calculate the interaction energy we resort to a well-known method ^Hj of 
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truncating the potential beyond a cut-off length l cu tof / 3> X a b in order to keep 
the computation time within a sensible limit. We make sure that our results 
do not depend on the size of lcutof f- Since we are interested in studying the 
effects of defect correlations in the dynamics, we restrict our simulation to low 
temperatures T/T* < 1. Here T* denotes the temperature where entropic 
corrections due to thermal fluctuations become relevant for pinned flux lines 
6 . The columnar pins and point defects are respectively modeled by uniform 
cylindrical or spherical potential wells of (uniform) strength U and radius bo . 
The cylindrical wells are oriented parallel to the c axis. We investigate the 
dynamics for three different distributions of defects, namely, (i) columnar pins 
either distributed randomly, or (ii) arranged in a square lattice in the xy plane, 
and (iii) point defects distributed randomly in the sample. Since in the presence 
of weak currents J = — Jy the flux lines locally move according to equilibrium 
dynamics, one may incorporate the effect of the force in the MC by introducing 
an additional work term — Jo ^ s ri ( s ) m ^ ne Haimltonian [51 1131 HI) . 

To initiate the simulation, flux lines of straight configurations parallel to the 
z axis are nucleated at random positions in the xy plane. We then let the system 
equilibrate. Once equilibrium has been reached we turn on the external current 
and let the system evolve. We have checked that the results in the steady state 
are independent of the initial configurations. At each trial a randomly chosen 
point on the line is updated according to a Metropolis algorithm 120] • The point 
can then move in a random direction a maximum distance of A < b /^/3 to 
guarantee interaction with every defect. In the simulation the values Zi are held 
fixed; we have checked that our results do not change qualitatively even if these 
positions Zi are allowed to fluctuate [To]. 

The drift velocity of the flux lines is proportional to the average velocity 
of the center of mass (CM) v cm = nT 1 Z)2= 1 ([Ricm('J") - Ricm(0)]/r), where 
Hi emir) — Ricm(O) is the distance traversed by the CM of the ith flux line in a 
time interval of length r; (• • •) and the overbar denote the averages over Monte 
Carlo steps (MCS) in the steady state and over different disorder realizations, 
respectively. The voltage drop measured in experiments is caused by the induced 
electric field E = B x v cm /c [21]. All length and energy scales are measured in 
units of bo and So, respectively. The average distance between the defects for 
a uniform random distribution and the lattice constant for the periodic array 
of columnar pins was taken as d = 15 bo- The parameters X a b, Cafc, £, U, and 
T were chosen to be 56 bo, 0.64 bo, 4e , 0.0075 eq, and 4 respectively in the 
simulation; these numbers are consistent with the experimental data for high- 
T c materials pQ. Therefore thermally induced bending and wandering of the 
flux lines are largely suppressed, and we can interpret our results in terms of 
low-temperature dynamics. In the simulations all sets of data were collected in 
the steady state, which was typically reached after t > 10 5 MCS and t > 10 6 
MCS for columnar defects and point pins, respectively. The size of the system 
ranged from 60 x 9 x 60 to 2000 x 9 x 60 in the simulations, and the data were 
averaged over 20 to 50 realizations of disorder. The value of r ranged from 30 
to 500 MCS in the simulations. 



4 



0.25 



V 



Figure 1: I-V characteristics for interacting flux lines in the presence of colum- 
nar defects with random (x) and periodic (o) spatial distribution. Note that J c 
is almost the same in both situations. The corresponding I-V plot for a single 
flux line with periodically arranged columnar defects (A) shows a higher J c (see 
text). In the non- interacting (single line) case with columnar pins J c is almost 
two times larger as compared to the interacting system. The I-V plot for a single 
flux line with random point defect distribution (*) is also given for comparison 
(see text). For currents J > J c the system with randomly distributed columnar 
pins displays a larger resistance (higher average vortex velocity) than the one 
with a periodic defect distribution. 

4 Results 



In the absence of disorder and any external drive, at low temperatures (in our 
simulations T < T* <C T c ) the flux lines arrange themselves in a triangular 
Abrikosov lattice |23 ■ in the presence of disorder the lattice structure is replaced 
by glass phases [5] , namely the so-called vortex glass in the case of uncorrelated 
point disorder, and the Bose glass in the presence of columnar defects 0. As 
an external current is applied the nearly localized vortices (since T < T* the 
flux lines rarely tunnel through the defects) in the glassy phases start moving 
once a critical driving current is reached. At T = 0, and in the absence of 
mutual vortex interactions, the flux lines are strictly bound to one or more pins 
in equilibrium. At the critical depinning current J = J c one then encounters a 
continuous non-equilibrium phase transition [221 1231 l2"l) . This zero-temperature 
phase transition becomes rounded to a crossover at T > 0, but a characteristic 
'critical' current separating the localized and moving vortex regimes can still 
be inferred at sufficiently low temperatures, as seen in the I-V characteristics 
displayed in Fig. ^ which were computed at T = 0.25 x 10~ 3 £ <C T*. 

It turns out that the critical current J c is actually higher for a system with 
columnar pins as compared to point disorder because of the increased defect 
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Figure 2: Log-log plot of the velocity noise power spectrum S(uj) vs. frequency 
u> for a system with columnar defects near the depinning threshold at J w J c . 
The solid line, given as a guide, has slope —2. 

correlations However the value of J c is the same for both random and peri- 
odic spatial distribution of the pins |16| . Repulsive interactions between the flux 
lines facilitate flux creep which is reflected in the derealization of the vortices 
at a considerably lower current than for the non-interacting case (Fig. ^| . The 
TV characteristics also reveal that the voltage at the same current is higher for 
the system with a random as compared to periodic defect distribution. Sim- 
ilar behavior is observed in the dynamics of the single flux line, and can be 
qualitatively explained by an effectively faster transit time in the random pin 
distribution compared to the periodic columnar defect array . We have also 
found that in the system with a single flux line in the presence of a periodic pin 
arrangement the I-V characteristics depends strongly on the orientation of J 
with respect to the lattice direction, for this determines the effective density of 
defects encountered by the moving vortex. A similar feature is likely to survive 
in the presence of interactions. 

The velocity (proportional to the induced voltage) noise is calculated by 
computing the velocity fluctuations about v cm . We evaluate the power spec- 
trum S(w) = v x (oj)v*(uj), where v x {ui) is the Fourier transform of the velocity 
fluctuation v x = n _1 X)"=i i Vix ~ witn v lx {t) = [X icm (t +r) -X jcro (r)]/r 

being the center-of-mass velocity of the ith flux line. 

Near J c we observe a power-law decay in the voltage noise spectrum: S(u>) ~ 
l/uj a , where an effective exponent a rs 2 is observed for almost a frequency 
decade in the system with periodically arranged columnar pins (Fig. 0). At 
low frequencies there are indications that the effective exponent may be smaller 
than 2. For a single line in the presence of random point defects and at T = 0, 
the value of a can be determined through a functional renormalization group 
calculation which gives a — 1.5 to one-loop order in three dimensions |2HI 12 lj . 
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Figure 3: (a) Plot of the velocity noise power spectrum S(u>) vs. uj/v cm for a 
system with columnar pins arranged in a periodic array for J — 0.05 (circles), 
J = 0.07 (squares), and J = 0.09 (triangles). The peaks occur at the same 
scaled frequency to cx v cm /d. Note that the width of the peak decreases as J 
is increased indicating the wider separation between the time scales r r and t v . 
(b) Plot of S(to) vs. uo for columnar pins distributed randomly in the sample 
for J = 0.05 (squares). Notice the presence of secondary peaks in the spectrum 
(see text). 

with corresponding mean-field value a = 2. The power law behavior observed in 
the simulation (Fig. [2J) can be interpreted as a remnant of the zero-temperature 
depinning transition. In experiments [Ij, such broad-band noise (BBN) features 
show a uj~ a decay with a w 2. 

For currents J > J c , three distinct time scales can appear in the dynamics: 
the residence time r r which is the duration a vortex line is trapped on the 
defects, t v , the time taken by the flux line to travel between the pins, and 
t w , the time taken by the moving vortex lattice in the event of reordering to 
traverse a distance of the lattice constant. In the regime where these time 
scales are well separated and r r < t v , as well as r r < r w , we can expect t v 
and t w to produce pronounced peaks in the velocity noise spectrum. We first 
summarize the results from the non-interacting case where the time scale t w is 
entirely absent (stick-slip motion). In Fig.^a) we plot S (u>) at different driving 
currents for a periodic distribution of columnar pins. The spectra reveal a single 
peak corresponding to the frequency Co cx v cm /d cx l/r„. Note that the width of 
the peak decreases as the external current J increases, which widens the time 
scale separation between r r and t v . In a random defect distribution, see Fig. 
Etb), there would be a distribution of the time scales t v around the average value 
d/v cm . Therefore there could be many secondary peaks present in the spectrum 
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Figure 4: Plot of the velocity noise spectrum S(lu) vs. Lu/v cm for a system of 
periodically arranged columnar pins at J = 0.07 with interacting flux lines. The 
primary peak corresponds to the rescaled frequency Q oc v cm /d. 



arising due to the distribution of t v in addition to the principle peak which 
corresponds to t t oc d/v cm - These secondary peaks would eventually diminish 
in intensity as the distribution of the time scales will narrow down at d/v cm 
with the results averaged over more defect configurations. Unlike for correlated 
defects, the time scales r r and t v are not well separated for point pins and the 
peaks turns out to be considerably suppressed as compared to the situation for 
columnar defects 

The presence of interactions between the flux lines may introduce the third 
time scale t w if the vortex lines align themselves into a lattice structure at 
large external currents. For example the predicted moving vortex lattice (in 
the presence of point defects) [7j would give rise to a distinct t w cx a m /v cm , 
where a m is the lattice constant of the moving flux line array. We have not yet 
detected any signature of the corresponding narrow-band noise in the presence 
of columnar defects with periodic spatial distribution, even when the average 
distance between the flux lines is much less than the periodicity of the pins. 
This implies absence of reordering in the observed range of external currents. 
We do find a peak in S(u>) corresponding to t v oc d/v cm (Fig. 0J) which is also 
observed in the non-interacting case. However we might expect reordering when 
the density of the flux lines is larger than the defect density because then the 
effect of defect correlations becomes masked. 

We have noticed a marked difference in S(oj) between columnar and point 
pins which results from the lack of spatial correlation in the z direction for the 
point disorder. Notice that this distinction can be observed only in a three- 
dimensional dynamical simulation. It is due to a generic difference between 
correlated and localized pinning potentials. The ensuing absence or presence of 
the narrow-band noise peaks in the voltage noise spectrum S(uj), and its specific 
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features, could perhaps be utilized as a signature to identify and characterize 
the type of disorder present and responsible for flux pinning in experimental 
superconducting samples as well. 

5 Conclusion 

In conclusion we have investigated the three-dimensional dynamics of interacting 
flux lines in the presence of correlated defects, either arranged in a periodic ar- 
ray or placed randomly in space. We find that the nature of the current- voltage 
characteristics depends on the defect correlation as well as the spatial distri- 
bution of the pinning centers. The voltage noise spectrum shows pronounced 
peaks in the case of correlated pins. The difference in the power spectra at 
J > J c between columnar and point disorder may serve as a novel diagnostic 
tool to characterize the pinning centers in superconducting samples. 
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